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The Hartree-Fock Perturbation Theory for Finite Dimensional Basis Sets. 
Part II. Construction and Application of the Variable Basis Sets in Calculations 
of Dipole and Quadrupole Polarizabilities 


1. WSTĘP 


Oddziaływania między dwoma układami wieloelektronowymi (atom, jon, molekuła) ` 
można rozpatrywać jako wynik oddziaływania rozkładu gęstości elektronowej jednego 
z nich z polem elektrycznym wytworzonym przez drugi [1]. Rozkład gęstości elektro- 
nowej reprezentowany być może przez trwałe, indukowane i fluktuacyjne momenty 
multipolowe [2]. Podobnie pole elektryczne pochodzić może od trwatych,indukowanych 
i fluktuacyjnych multipoli [1—3]. Przy założeniu dużej odległości R między układami 
mamy stąd trzy typy oddziaływań: elektrostatyczne, wynikające z oddziaływania trwa- 
łych momentów multipolowych obu układów, indukcyjne, będące wynikiem oddziały- 
wania trwałych momentów jednego układu z momentami indukowanymi w drugim ukła- 
dzie, oraz dyspersyjne, pochodzące z oddziaływania momentów fluktuacyjnych obu — 
układów [1—3]. 

Wielkość trwałych momentów multipolowych determinowana jest wyłącznie przez 
strukturę rozkładu gęstości elektronowej danego układu. O wielkości momentów indu- 
kowanych i fluktuacyjnych decyduje ponadto charakter i natężenie zewnętrznego pola 
elektrycznego. W słabym polu elektrycznym wielkość indukowanego momentu multipo- 
lowego jest wprost proporcjonalna do natężenia tego pola, a współczynnik proporcjonal- 

* Praca była częściowo finansowana prżez Instytut Niskich Temperatur i Badań Strukturalnych 
PAN w ramach tematu nr MR-1.9.4.3/2. 
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ności, wynikający ze struktury rozkładu ładunku elektrycznego danego układu, nosi 
nazwę polaryzowalności multipolowej [1, 2, 4, 5]. Za pomocą tych wielkości oraz trwa- 
‘tych momentów multipolowych można wyrazić dokładnie energię indukcyjnego oddzia- 
ływania między układami [1—3]. Charakter energii dyspersyjnej jest bardziej złożony 
[1, 2,6]. Tym niemniej w przybliżeniu Unsólda [7] można ją wyrazić poprzez polaryzo- 
walności multipolowe oddziałujących układów. Tak więc pełna energia odziaływania, 
będąca sumą energii elektrostatycznej, indukcyjnej i dyspersyjnej, może być wyrażona 
poprzez trwałe multipole i polaryzowalności multipolowe izolowanych układów. 

Trwałe momenty multipolowe oraz połaryzowalności multipolowe są charakterystycz- 
ną własnością danego, izolowanego układu wieloelektronowego. Reprezentują one odpo- 
wiednio rozkład gęstości elektronowej układu i jego podatność na deformację pod wpły- 
wem zewnętrznego pola. 

Polaryzowalności multipolowe determinują wielkość indukcynej i dyspersyjnej energii 
oditziatywan, które można przedstawić jako nieskończone sumy postaci [6, 12]: 


Cn 


DS n=4 R” dysp n= RV 


Współczynniki A zawierają iloczyny trwałych momentów i polaryzowalności oddzia- 
łujących układów, a współczynniki Co można w przybliżeniu przedstawić jako iloczyny 
ich polaryzowalności multipolowych. Z uwagi na założoną dużą wartość R istotny wkład 
- do energii indukcyjnej i dyspersyjnej mają pierwsze człony obu szeregów, determinowane 
przez najniższe polaryzowalności multipolowe, to jest dipolową i kwadrupołową, Polary- 
zowalności te występują w energii indukcyjnej, począwszy od członów odpowiednio rzędu 
R”* iR”5,awenergii dyspersyjnej R"* i R“? [6, 12]. 

Doświadczalne metody wyznaczania polaryzowalności dipołowych oparte są na bada- 
niach refrakcji optycznej lub polaryzacji elektrycznej [3] zjawiska Kerra [8,9], oraz po- 
- miarze anizotropowego rozpraszania światła [10, 11]. Z uwagi na subtelność efektów bę- 
dących podstawą pomiaru uzyskiwane wyniki są nie zawsze zadowalającej jakości [1, 2, 
12]. Zdecydowanie mniejsze są obecnie możliwości doświadczalnego wyznaczenia polary- 
zowalności kwadrupolowych [12]. W tej sytuacji dużego znaczenia nabierają teoretyczne 
metody obliczania tych wielkości, oparte w większości na zastosowaniu rachunku zabu- 
rzeń w przybliżeniu Hartree-Focka [13—15]. Praktyczne wykorzystanie metody CHF 
wymaga jednak zastąpienia dokładnych rozwiązań Hartree-Focka przez przybliżone 
rozwiązania SCF uzyskane w skończenie wymiarowej bazie funkcyjnej [15]. Pojawia się 
wówczas problem zbieżności otrzymywanych wyników do wartości dokładnych mętody 
CHF [15, 16]. W tradycyjnych obliczeniach metodą CHF w ustalonych bazach funkcyj- 
nych zbieżność ta jest bardzo wolna, co pociąga za sobą konieczność stosowania baz 
o dużych wymiarach [17]. To z kolei powoduje wzrost trudności numerycznych [15], 
które mogą ulec istotnej redukcji, jeżeli w obliczeniach zastosowane zostaną bazy funkcji 
nie ustalonych, lecz bezpośrednio zależnych od zaburzenia [18—24]. Możliwość poprawy 
zbieżności w metodzie CHF przez zastosowanie takich zmiennych baz funkcyjnych była 
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omawiana w części I*. W części II omówiony zostanie sposób konstrukcji zmiennych baz 
‘funkcyjnych odpowiednich dla obliczeń polaryzowalności dipolowych i kwadrupolowych. 


_ 2. POLARYZOWALNOSCI MULTIPOLOWE 


Polaryzowalności multipolowe są wielkościami tensorowymi rzędu 2/ [1, 2]. Dla ko- 
lejnych wartości / = 1, 2, ... otrzymamy stąd tensory : drugiego rzędu polaryzowalności 
dipolowej a; , czwartego rzędu polaryzowalności kwadrupolowej a, . itd. Omawiane tutaj 
wielkości wygodnie jest zdefiniować, wprowadzając zewnętrzne pole elektryczne o stałym 
(1 — 1) gradiencie natężenia (tensor /-tego rzędu). Energia danego układu w tym ze- 
wnętrznym polu elektrycznym może być zapisana jako** 


1 . 
EV) =E% =g Mu Eé CC KE 


gdzie E (0) oznacza energię układu izolowanego, a Hr jest jego momentem multipolowym 
l-tego rzędu. Wszystkie występujące w tym równaniu tensory (py, a, V,) można wyrazić 
bądź we współrzędnych kartezjańskich, bądź we współrzędnych sferycznych [1, 25]. 
W tej pracy zastosowano układ współrzędnych sferycznych. 

Powyższe równanie pozwala zdefiniować składowe momentów i polaryzowalności 
multipolowych odpowiednio poprzez pierwsze i drugie pochodne energii E (V; ) względem 
składowych gradientu pola V 


dE(V)) 3 E(V)) 
RE lest, VEER d. th = e 
| Fm |V=0 TELE 
gdzie m,n = —1, + l. Tak zdefiniowane wielkości ki a, są identyczne z momentami 


i polaryzowalnościami multipolowymi wprowadzonymi przez Muldera [25]. 
Ostatnie równania ukazują związek omawianych wielkości z rachunkiem zaburzeń. 
Mamy mianowicie: 
kim = EU, >  Gd;mn = RRE 
Wielkości Æ} (1) oraz E} kel n 54 Odpowiednio energiami pierwszego i drugiego rzędu 
w rachunku zaburzeń: ; 


* K. Woliński: Ann. Univ Mariae Curie-Skłodowska, Lublin, Sectio AA, 34 (1979). 
` ** We wzorze występuje mnożenie tensorów odpowiedniego rzędu. 
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Funkcja y © olka stan układu réezaburzonego, a funkcje RI i v(t) są funkcja- 
mi pierwszego rzędu względem zaburzenia odpowiednio —V/7., fi]; m oraz —Vj-n fin, 
Operator jij: jest operatorem m-tej składowej momentu multipolowego układu. Jak 
widać, obliczenie pozadiagonalnych składowych polaryzowalności multipolowych 
(n +m) wymaga stosowania rachunku podwójnych zaburzeń [22, 26]. 

_ Dla układów o symetrii sferycznej, z uwagi na izotropię oi. an = Bang, polary- 
zeswalność multipolowa jest określona przez jedną składową (np. a7;90) równą jedno- 
cześnie średniej PCA multipolowej a; 


1 ` mes 
el EECHER E POK 
Wielkość tę można wyrazić poprzez energię drugiego rzędu Ef2), = EI) 
a, = AEN l 
w rachunku ii e z pojedynczą perturbacją reprezentowaną przez 
H® = — Vue Bug. 


Operator składowej momentu multipolowego /i 1:0 jest sumą operatotów jednoelektro- 
nowych postaci: 


SI? Ste? (cos 9) . 


PEATA RAEC do przypadku EE eh (l = 1) i kwa- 
drupolowych (1 = 2), otrzymamy: 


a = 04,00 = — 2542) eg = Q2;00 = —2E$2) , 


gdzie dla obu przypadków jednoelektronowe operatory zaburzenia mają odpowiednio 
postać 


rcosQ = z 


Q1,0 | 
am Zen cov@ — 1) = Ca —x?—y?), 


W nastepnym punkcie tej pracy omówiony zostanie sposób konstrukcji baz funkcyj- 
nych bezpośrednio zależnych od zaburzenia —V 02 1:0 da 1 = 1 i 7 = 2. Bazy takie 
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były zastosowane do obliczeń dipolowych polaryzowalności atomów i molekuł [19—22] 
- oraz ostatnio do obliczeń polaryzowalności kwadrupolowych [23, 24]. 


` 3. ZMIENNE BAZY FUNKCYJNE 


Koncepcja użycia zmiennych baz funkcyjnych w obliczeniach elektrycznych włas- 
ności atomów wprowadzona została przez Moccia [18]. Zaproponował on sposób 
* konstrukcji zmiennych, baz funkcyjnych odpowiednich dla obliczeń polaryzowalności 
dipolowych. Mając bazę funkcji dla problemu niezaburzonego dea ]. utworzył on bazę 
funkcji zależnych od zaburzenia — V1 ; o Q1;0, mnożąc każdą funkcję fi przez wspólny 
czynnik wykładniczy: 


Mo) = AE exp (Vigo a. 


Dodatkowy parametr k wyznaczany był z warunku minimum energii drugiego rzędu. 

Wynik uzyskany przez Moccia dla dipolowej połaryzowalności atomu He —1.1 au. 
(z k = 0.68) był jednak dość daleki od dokładnej wartości Hartree-Focka —1.32 au. [27]. 
Powodem tego niepowodzenia były prawdopodobnie niedpowiednie własności funkcji 
JE: el 

Zaproponowany przez Moccia sposób konstrukcji zmiennych baz funkcyjnych nie 
ma prostego uzasadnienia fizycznego. Sposób konstrukcji, który zostanie przedstawiony 
niżej, oparty jest natomiast na rozważaniu wpływu zaburzeń typu —V;;9Q7;0 (Z 1 = 
= 1, 2) na trójwymiarowy oscylator harmoniczny. 

Rozwiązaniem zagadnienia własnego dla niezaburzonego_ PM harmonicznego: 


[- ŻY? +20 ]g = "e 


są, jak wiadomo, iloczyny funkcji Gaussa i wielomianów Hermite'a, które można zapisać 
w postaci; 


; & = Segel" Sy (Vay) * gz(Ziaz), 


gdzie a = ay = ay = az jest wykładnikiem funkcji g. Właśnie fakt, że rozwiązania dla 
oscylatora harmonicznego są związane z funkcjami Gaussa, spowodował zainteresowanie 
tym modelem fizycznym. Funkcje gaussowskie są bowiem często stosowane jako bazy 
w obliczeniach kwantowo-chemicznych, 

Jeżeli oscylator poddany zostanie zaburzeniu postaci —V7;9 Q7:9 Z 1 = 1, 2,to rozwią: 
zaniem równania: 


b” ; V? + 2a°r? — Viso Qizo] s(V1;o) = e(Vi;0) &(Vi;0) 


sa odpowiednio zmody fikowane funkcje g. Dla / = 1 można je zapisać jako: 
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g(V1;0) = &x (idx) * 8y Wray) * Se (el: gel, 
i b 

Gat "ue 

Dla przypadku / = 2 otrzymamy natomiast: 


gdzie z =z — 


DCH Ee E , 
ZE: 
4a SE Va; E? 
Tak więc zaburzenie o charakterze dipolowym (/ = 1) powoduje w przypadku funkcji 
g ogólnej postaci g(r;R,a) przesunięcie jej centrum R = [X, Y, Zk 


-F A; ` 
R Word i R (V1;0), 


gdzie ax = ay = a(1 + [> „da z Kz; 2 ia, = a(1 — 


gdzie R' = R + AR z AR; = 5; SS Haf = x,7,2) [19]. Nie ulega przy tym zmianie 
a 


sferyczna symetria wykładniczej części funkcji 


ealt- R)? zaburzenie eot- R? 
=] p b 
Odwrotnie jest w przypadku zaburzenia o charakterze kwadrupolowym (/ = 2): poło- 
żenie centrum nie ulega przesunięciu, zmienia się natomiast wykładnik: 
Vaio Q330 


a'(V2;0) ` (2) 
gdzie: 
a; = = a[1 ek (1 WEE 7 V2; Ei 


Powoduje to zmianę symetrii wykładniczej części „łe ze sferycznej na eliptyczną: 


e-a(r-R)? _ zaburzenie g 2x (x - X)* -ay(y- Y) -az(z-Z)? ` 


r="2 


Powyższe rozważania pozwoliły ustalić fizycznie uzasadnioną postać zależności funkcji 
Gaussa od zaburzenia —V;.g Q7;9 dla / = 1 il = 2. Mając daną bazę funkcji Gaussa dla 
problemu niezaburzonego można w oparciu o te zależności wyprowadzić zbiór funkcji 
zależnych od V;-9 oraz zbiór funkcji zależnych od V2.9. Funkcje te opisują stan zbioru 
niezależnych oscylatorów w zewnętrznym polu elektrycznym o stałym odpowiednio 
natężeniu i gradiencie natężenia. Ponieważ dany układ wieloelektronowy nie może być 
dokładnie reprezentowany przez zbiór niezależnych oscylatorów, to aby uwzględnić 
istniejące różnice wprowadza się dodatkowy parametr skalujący wpływ zaburzenia V;.g. | 
Parametr ten wyznaczany będzie z warunku minimum odpowiedniej energii drugiego ` 
rzędu. Dla przypadku dipolowego będzie to parametr k, [19]: 
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= 3 
3 6) 

a dla przypadku kwadrupolowego parametr ką: 
e a =.a[1+(1 WEE Lat. (4) 


Ostatecznie więc zmienne bazy funkcji ARGIA (GTO = Gaussian Type Orbitals), 
* zależnych od zaburzenia postaci —V;;9 Q;;0,wyprowadza się z początkowej bazy funkcyj- 
nej dla problemu niezaburzonego: f 


sit) = sin R ai | RAUT oki) = ALE RV; 0,K1),a] 
NZ ;0» Ka) = ALE R,a "Wa; 0>%2)] 


(5) 


gdzie R' ia’ określone są wzorąmi (3) i (4). 
Koncepcję uzmienniania baz GTO rozszerzono na bazy funkcji typu Slatera (STO = 
Slater Type Orbitals). W tym celu wykorzystano znaną transformację całkową [28], 
pozwalającą wyrazić wykładniczą część orbitalu STO przez wykładniczą część orbi- 

talu GTO: i 
RA 


,  exp(—dr) = Pai ax 2 exp (- 


b? 
) exp (—ar*) da, 
a > 


Podstawiając w miejsce wykładniczej części orbitalu GTO uzyskane poprzednio wyra- 
żenie, opisujące jej zależność od wielkości zaburzenia V;;9, otrzymamy równanie opisują- 
ce zmianę wykładniczej części orbitalu STO pod wpływem tego zaburzenia. Mając takie - 
wyrażenie, można już otrzymać zależne od zaburzenia —V 7.9 Qj;9 bazy STO: 


z Size, kı) = s(r;R,b, Ho, ki) 
so) = sj (GR,b) | 
si(Vi;0 k2) = s;GGR,b, Va ;0>K2) . 


Przedstawione wyżej zmienne bazy GTO i STO mają szereg ciekawych własności. 
Po pierwsze, wymiar zmiennej bazy jest taki sam, jak wymiar bazy początkowej (dla 
problemu niezaburzonego), ponieważ każda zmienna funkcja wyprowadzana jest z jednej 
funkcji początkowej. Po drugie, używając w obliczeniach energii drugiego rzędu, np. 
m-wymiarowej zmiennej bazy, korzystamy w rzeczywistości ze znacznie większej liczby , 
funkcji. Obliczając bowiem E f2), zgodnie z formułą podaną w Części I, posługujemy się 
funkcjami bazy w pierwszym i drugim rzędzie, które to funkcje są liniowymi kombinacja- 
mi, nie występujących w bazie początkowej, orbitali o wyższej symetrii. Na przykład, 
jeśli początkowa baza GTO zawierała orbital typu s, to wyprowadzony z niego orbital 
zależny od V;.9 będzie miał w pierwszym i drugim rzędzie postać: 
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BEE dla [= 1 [22] 
BP Hie dl + e, däi 


oraz _ 


> dla/=2 
DI = k z kå [cy gf) + cą dz?) + Ca 21. 
Ga, orbitale sf°), pf), df), gf) mają symetrię typu odpowiednjo s, p, d, 8g, a wykład- 
niki w nich występujące są takie same jak w orbitalu ff 0) bazy początkowej. Tak więc, 
korzystając :z funkcji bazy w pierwszym i drugim AS rozszerzamy w rzeczywistości 
bazę funkcyjną o wyższe orbitale, których parametry są dokładnie znane i nie wymagają 
dodatkowych optymalizacji. Dodać przy tym należy, że zarówno typ tych ,,dodatko- 
wych” orbitali, jak również ich wkład do funkcji pierwszego i drugiego rzędu jest wew- 
nętrznie kontrolowany przez samą metodę. Po trzecie, dzięki przedstawionemu sposo- 
„ bowi wprowadzania parametrów skalujących k; i ką odpowiednie energie drugiegi rzędu 

_ są parabolami względem nich. Pozwala to wyznaczyć w bardzo prosty sposób optymalne 
wartości ky i k3. ` 

Na zakończenie warto podkreślić, że największy wkład do energii drugiego rzędu be: 

dą miały orbitale z małymi wykładnikami. Zgodnie bowiem z równaniami (3) i (4) takie 
właśnie orbitale doznają największych zmian pod wpływem zewnętrznego zaburzenia. 
Jest to zrozumiałe pod względem fizycznym, ponieważ orbitale z małymi wyktadnikami 
dają główny wkład do gęstości elektronowej dla większych odległości od jądra. Gęstość 
elektronowa w tym obszarze decyduje z kolei o polaryzowalności układu, a tym R 
o wartości odpowiedniej energii drugiego rzędu. 


4. OBLICZENIA POLARYZOWALNOŚCI DIPOLOWYCH I KWADRUPOLOWYCH 


Opisane w poprzednim punkcie zmienne bazy GTO i STO zastosowane zostały do 
obliczeń polaryzowalności dipolowych [19—22] oraz kwadrupolowych. Odpowiednie 
energie drugiego rzędu obliczane były metodą CHF zgodnie z równaniem na £(?), poda- 
nym w Części I. Optymalne wartości parametrów E, i ką wyznaczane były analitycznie 
z warunku minimum paraboli E{?)(k,) i £§?)(k2). Zmienne bazy. GTO i STO wyprowa- 
dzane były z początkowych, optymalizowanych względem E (0), baz GTO vanDuijne- - 
veldta [29] oraz baz STO Clementiego [30]. g 

Pierwsze testowe obliczenia przeprowadzone zostały da atomu wodoru. Układ ten 
` jest bardzo wygodny do prowadzenia badań testowych, ponieważ znane są w tym wy- 
padku dokładne rozwiązania zarówno dla problemu niezaburzonego, jak również zabu- 
rzonego z perturbacjami typu —F;:9Q7:0 [4]. Wee m rozwiązaniem dla układu 
niezaburzonego w stanie podstawowym jest funkcja f/©) = N exp (— 71). Dla zaburzeń 
omawianego tutaj typu dokładne funkcje pierwszego rzędu mają postać: 
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1 
Kx = Guer FOI, Fx = Zd + ZO. 


Zmienne orbitale STO wyprowadzone z funkcji fo zgodnie z opisaną poprzednio 
metodą dają w pierwszym rzędzie funkcje postaci: 


ję = ky Qio (l+ NF, £0) = Zb Q;0(1 + NSO. 
Optymalne wartości parametrów k, i ką wyznaczone z minimum paraboli: 
EG) = 88 k? —28k,, E$) = 13.6875 k — 20.25 k 


są równe odpowiednio 0.636363 oraz 0.739726. Odpowiadają im minima E$2) = 
= — 2.2273 [20] oraz-E{?) = — 7.489726, które pozostają w bardzo dobrej zgodności 
z dokładnymi wartościami Ef?) „ = — 2.25 i E$?) „ = — 7.5. Równie dokładne rezulta- 
ty otrzymano dla atomów i jonów wieloelektronowych. Przedstawiono je w tab. 1, w któ- 
rej obok wyników otrzymanych metodą CHF w zmiennych bazach funkcyjnych zamiesz- 
czono również rezultaty uzyskane przez innych autorów w ustalonych bazach o dużych 
wymiarach. 

„Dokładność wyników otrzymanych omawianą metodą w istotny sposób zależy od 
? wymiaru i jakości bazy funkcyjnej dla układu niezaburzonego, z której wyprowadzana 

jest baza zależna od zaburzenia. Obrazuje to tab. 2. Zamieszczono w niej wyniki dla dipo- 
lowej i kwadrupolowej polaryzowalności atomu He, uzyskane w różnych bazach GTO 
i STO. Jak widać, podstawą uzyskania dokładnych wartości polaryzowalności jest użycie 
dobrych (ze względu na E(%)) baz funkcyjnych dla układu niezaburzonego. Tabela 2 po- 


r Tab. 1. Dipolowa i kwadrupolowa polaryzowalność atomów i jonów 
e (wszystkie wielkości w jednostkach atomowych) 


Polaryzowalność 
Atom 7 
jon dipolowa kwadrupolowa 
GTO? ` snob ` eg - Grof sto! CHF 

He 1.3214 1.3193 ` 1.3227 2.323 2.326 2.330 
Lit 0.1889 0.1896. 0.112 0.113 
Be*? 0.0514 0.0519 0.015 0.016 
Bt 0.0192 0.0196 0.003 0.004 
Gt wz 0.0089 0.0089 0.001 0.001 
Be 43.344 45.6° 340.76 339.55 347.02 
- B Ref. 22. 

b A.J. Sadlej: Acta Phys. Pol. A53, 297 (1978). 

© S.T. Epstein, R. E. Johnson: J. Chem. Phys. 47, 2275 (1967). 

d P. W. Langhoff, M. Karplus, R.P. Hurst:J. Chem. Phys. 44, 505 (1966). 

ż R. R. Teachout, R. T. Pack: Atomic Data 3, 195. 


Wyniki własne. 


` 
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Tab. 2. Zbieżność w metodzie CHF w zmiennych bazach 
GTO i STO. Dipolowa i kwadrupolowa polaryzowalność 
atomu He w jednostkach atomowych 


Lee r AER, D ab A 


2 2.747066 0.8834 0.701 u.936 1.092 
3 2.835680 1.0866 0.634 1.424 0.940 
5 2.859896 1.2720 0.571 2.050 0.800 
8 2.861625 1.3177 0.544 2.290 0.733 
10 2.861673 1.3214 0.539 2.316 0.723 
T 
1 
2 
4 
5 


2.8476563 1.1226 0.6504 3.642 2.077 

2.8616700 1.3210 0.5404 2.330 0.726 
2.8616785 .1.3193 0.5328 2.256 0.676 i 
2.8616799 1.3212 0.5360 2.334 0.707 


a Ref. 22. 
Wyniki własne. 


kazuje również, że zmienne bazy STO są bardziej efektywne niż bazy GTO, głównie jed- 
nak w przypadku polaryzowalności dipolowych, gdzie wynik uzyskany w bazie 2 STO 
i bazie 10 GTO jest praktycznie taki sam. W obliczeniach polaryzowalności kwadrupolo- 
wej przewaga baz STO nad bazami GTO jest mniejsza. 

Stan podstawowy izolowanych atomów i jonów (H, He—C**, Be) może być opisany 
w bazach zawierających wyłącznie funkcje typu s. Metoda CHF w ustalonych bazach 
funkcji typu s daje dipolowe i kwadrupolowe polaryzowalności równe tożsamościowo 
zeru, co wynika z symetrii odpowiednich operatorów zaburzenia. Natomiast użycie zmien- 
nych baz wyprowadzonych z funkcji typu s pozwala uzyskać wyniki o bardzo dobrej dò- 

kładności, przy czym konieczna jest tylko prosta optymalizacja parametru ky. 

Na zakończenie warto podkreślić, że zastosowania przedstawionej metody nie ograni- 
czają się do wyznaczenia polaryzowalności dipolowych i kwadrupolowych. Możliwe jest 
bowiem obliczenie wszystkich tych wielkości, które są związane z kolejnymi, wyższymi 
poprawkami energetycznymi w rachunku zaburzeń z perturbacją typu — Hr. m Q1;m dla 
l= 1il= 2. Są to, jak wiadomo, dipolowe i kwadrupolowe hiperpolaryzowalności róż- 
nych rzędów. Ponadto metoda nadaje sde do obliczania dipolowych i kwadrupolowych ` 
czynników ekranowania jąder. Te ostatnie obliczenia są obecnie prowadzone. Pierwsze 
otrzymane wyniki, szczególnie dla czynników kwadrupolowych, są wysoce zadowalające. 


Pragnę serdecznie podziękować Panu Doc. Drowi hab. Andrzejowi J. Sadlejowi za ree PRE 
okazaną przy pisaniu tej pracy. 
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PE3OME 


B- manno pa6oTe OÓcyxXJIeHo METOM IIOCTpo€HHA MepemMeHHEIX Ga3HCOB QyHKIHH, KOTOPKIE ` 
HeliocpecTBeHHO 3ABHCAT OT HaIIpAxeHHA, A TAKXKE OT rpapMeHTa HaIpAXKCHHA ANEKTPHYECKOTO ` 
nons. IlepemeHHkie QyEKIAt GBDIM NOJIYgEHBI HA OCHOBE AHAJIM3A COOTBETCTBCHHBIX pELIEHKK ZLA: 
TapMOHHYECKOTO OCUHJDIATOpa. JaTeM, DepemeHHbie Oa3HCHI dn HKUHH Dn 1 pHJIOXKCHBI DO SR 
NEHHA KBANDYNANBHBIX UOJAJA3YEMOCTK ATOMOB. OGHapy;KeHO Dome RE JITEX 
nepeMeHHbIX 6a3HCOB. 
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SUMMARY 


The method of constructing the GTO and STO bases which depend explicitly on the electric field 
strength and the electric field gradient strength is discussed. The appropriate variable basis sets are 
derived from the solution for the perturbed harmonic oscillator model. They have been applied to the 
calculation of atomic quadrupole polarizabilities in the framework of CHF method. Similar results 
obtained previously for the dipole polarizabilities are also analysed. It is pointed out that the proposed 
variable basis sets are very efficient in polarizability calculations. Moreover, other possibilities of 
application of the present approach are also discussed. 


